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ABSTRACT 

By correcting the redshift-dependent distances for peculiar motions through some pecuhar velocity 
field models, we recover the true distances of an extensive, all-sky sample of nearby galaxies 6400 
galaxies with recession velocities cz < 5500 km s^^ ), which is complete up to the apparent limiting 
blue magnitude J5 = 14 mag. Relying on available catalogs of galaxy groups, we treat ~2700 objects as 
members of galaxy groups and the remaining objects as field galaxies. 

We invert the derived redshift-distance relations to estimate distances for field galaxies and groups and 
we overcome the ambiguity inherent to the triple- valued zones by using Tully-Fisher relations calibrated 
on suitably defined samples of galaxies having distances predicted by peculiar velocity models. 

We use two independent approaches to modeling the peculiar velocity field: i) a cluster dipole re- 
construction scheme that we modify with the inclusion of a local model of Virgocentric-infall; ii) a 
multi-attractor model fitted to the Mark II and Mark III catalogs of galaxy peculiar velocities. In the 
multi-attractor model we assume that the velocity field is generated by a few prominent gravitational 
sources (Virgo cluster, "Great Attractor", Perseus-Pisces Supercluster, and Shapley supercluster). 

We discuss differences in the results coming out from different velocity models and from different Mark 
II and Mark III data subsets. In particular, according to Mark III data the Great Attractor appears 
to have a smaller influence on local dynamics than previously believed, whereas the Perseus-Pisces and 
Shapley superclusters acquire a specific dynamical role. Remarkably, the Shapley structure, which is 
found to account for nearly half the peculiar motion of the Local Group, is placed by Mark III data 
closer to the zone of avoidance with respect to its optical position. On the other hand, the modified 
cluster dipole model is characterized by relatively small flows towards the Great Attractor and the 
Shapley supercluster, together with a large Virgocentric infall. 

Our multi-attractor model based on Mark III data favors a cosmological density parameter flo ~0.5 
(irrespective of a biasing factor of order unity). 

The use of different peculiar velocity field models allows us to check to what extent differences in 
current views on cosmic flows affect the recovering of galaxy distances. We find that differences among 
distance estimates are less pronounced in the ~ 2000 — 4000 km distance range than at larger or 
smaller distances. In the last regions these differences have a serious impact on the 3D maps of the 
galaxy distribution and on the local galaxy density — on small scales (< 1 Mpc) — , which is a crucial 
parameter being used in statistical studies of environmental effects on the properties of nearby galaxies. 

Subject headings: galaxies: distances and redshifts — galaxies: clusters: general — cosmology: large- 
scale structure of universe 



1 INTRODUCTION 

The determination of galaxy distances is one of the 
most fundamental problem in extragalactic astronomy. In 
particular, galaxy distances allow us to map the three- 
dimensional (3D) distribution of galaxies and hence to 
evaluate the local galaxy density, which is an important 
characteristics of the galaxy environment. 

In many previous studies which rely only on the two- 
dimensional (2D) projection of galaxies on the sky (see, 
e.g., Lahav (1987) for all-sky 2D optical maps of bright 
galaxies), a local galaxy density parameter is defined in 
terms of galaxy counts over a sky area which are statisti- 
cally corrected for projection effects. With the advent of 
large surveys of galaxy redshifts and well-selected galaxy 



catalogs, it has become possible to map the spatial distri- 
bution of galaxies and hence to construct a 3D definition 
of the local galaxy density. 

Optical galaxy samples are more suitable for mapping 
the galaxy density field on relatively small scales than 
IRAS-selected galaxy samples, which have been frequently 
used as tracers of the density field on large scales up to 
a recession velocity cz ~ 20000 km s~^ (e.g., Rowan- 
Robinson et al., 1990; Kaiser et al., 1991; Strauss et al., 
1992a,b; Fisher et al., 1995). The latter samples have the 
advantage that they are homogeneously selected and pen- 
etrate to lower galactic latitudes because IRAS fluxes are 
much less impeded by Galactic extinction. But the IRAS 
galaxy samples have the serious drawback that they do not 
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include the early-type galaxies (because of their little dust 
content and star formation), which usually lie in the dens- 
est regions of clusters and superclusters. Furthermore, the 
density field traced by IRAS galaxies is noisy nearby, be- 
cause of the sparseness of the samples and because IRAS 
fluxes are much less linked with galaxy mass than optical 
fluxes. Therefore, in the following we shall focus on all-sky 
optical galaxy samples. 

The magnitude-limited Revised Shapley-Ames (RSA) 
Catalog of Bright Galaxies (Sandagc & Tamman, 1981), 
which contains the first large compilation of redshift data 
for bright galaxies across the entire sky, was used by Yahil, 
Sandagc & Tamman (1980) to delineate the 3D-dcnsity 
field of galaxies in the Local Supcrcluster (LS). The struc- 
tures of the LS region have been comprehensively de- 
scribed and named by TuUy & Fisher (1987). The former 
author intended to include all nearby galaxies with sys- 
temic recession velocities cz < 3000 km (2367 objects) 
in his Nearby Galaxies Catalog (Tully, 1988a. hereafter 
NBG), which is a combination of the RSA catalog and a 
diameter-limited sample of late-type and fainter galaxies 
found in an all-sky HI survey made principally by Fisher 
& Tully (1981) and Rcif ct al. (1982). In the NBG cat- 
alog, which is complete up to the corrected blue total 
magnitude Bt ~ 12 mag (although it extends to fainter 
magnitudes), the distances of all non-cluster galaxies have 
been estimated on the basis of their redshifts, with an as- 
sumed Hubble constant Hq = 75 kms~^Mpc^^ , whereas 
the galaxy members of systems with relatively high veloc- 
ity dispersion have been given a distance consistent with 
the mean redshift of the system. The important fact that 
redshifts are not equivalent to distances because galaxies 
have peculiar motions (i.e. deviations from the pure Hub- 
ble flow) is simply taken into account in the NBG catalog 
through the Virgocentric retardation model described by 
Tully & Shaya (1984), in which the authors assume that 
the Milky Way is retarded by 300 km s~^ from the pure 
Hubble flow by the mass of the Virgo cluster. Defining 
the local galaxy density on the basis of his NBG catalog, 
Tully (1988b) incorporated corrections for the catalog in- 
completion at large distances. 

Later, similar quantifications of the local galaxy density, 
based on the NBG catalog data, have been exploited in 
statistical analyses of environmental effects on some prop- 
erties of LS galaxies such as the frequency of bars (Giuricin 
et al., 1993), the classification of spiral arms (Giuricin et 
al., 1994), the frequency of active galactic nuclei (LINER 
and Seyfert objects) selected from optical spectroscopic 
surveys (Monaco et al., 1994), and the bulge-to-disk light 
ratio (Giuricin et al., 1995). 

In an effort of going beyond the LS, Hudson (1993a, b; 
1994a, b) assembled an extensive galaxy sample from a 
merging of the diameter-limited northern UGC catalog 
(Nilson, 1973) and the diameter- limited southern ESO 
catalog (Lauberts, 1982; Lauberts & Valentijn, 1989). He 
took into account the fairly large incompleteness in red- 
shift of his sample as a function of angular diameter and 
position on the sky, by applying statistical corrections, 



which allowed him to reconstruct the density field of opti- 
cal galaxies to a depth of cz = 8000 km s~^ . 

The Lyon-Meudon Extragalactic Database (LEDA), 
which collects and homogenizes several data for all the 
galaxies of the main optical catalogues, such as the cata- 
logs UGC (Nilson, 1973), ESO (Lauberts, 1982, Lauberts 
& Valentijn, 1989), CGCG (Zwicky et al., 1961-1968), and 
ESGC (Corwin & Skiff, 1998, in preparation), allows the 
extraction of an all-sky galaxy sample which has properties 
of completeness and is deep enough to probe the galaxy 
field significantly beyond the LS. In particular, Garcia et 
al. (1993) extracted a magnitude-limited sample which 
covers all the sky and contains 6392 galaxies having re- 
cession velocities cz < 5500 km s^^ . Although differ- 
ent optical catalogs are characterized by different limits of 
completeness in apparent magnitude or angular diameter, 
the above-mentioned galaxy sample was found to be sub- 
stantially complete up to its limiting corrected total blue 
magnitude Bt = 14 mag (for galactic latitudes |fe| > 20°) 
(Garcia ct al., 1993). The authors tabulated several pa- 
rameters for each galaxy, such as the morphological type, 
the maximum velocity rotation deduced from the depro- 
jectcd 21-cm hydrogen line width, the de Vaucouleurs' lu- 
minosity index, the distance modulus, the corrected an- 
gular sizes and total blue magnitudes. The authors de- 
rived the last two quantities by transforming the original 
raw data to the standard system of the RC3 catalog (de 
Vaucouleurs et al., 1991) and by applying corrections for 
Galactic extinction, internal extinction, and K-dimming. 
Only for the few galaxies which lack magnitude measures, 
the values of Bt are rough estimates derived from their 
diameters through a mean magnitude-diameter relation. 
The tabulated distance modulus was a weighted mean 
derived from the combination of redshift-distances (with 
Ho = 75 km s~^Mpc~^) and redshift-independent dis- 
tances obtained from distance indicators (DIs), i.e. from 
a TuUy-Fisher (TF) relation (which relates the absolute 
magnitude Mb to the logarithmic maximum rotation ve- 
locity logV^n) and a luminosity-luminosity index relation. 
In particular, the authors used the former kind of distances 
alone for the farthest galaxies (with cz > 1500 km s~^ ) 
and the latter kind alone for the few nearest objects (with 
cz < 500 km s^i ). 

Group assignments for the galaxies of this sample have 
been provided by Garcia (1983). She employed two 3D 
methods of group identification, the percolation or friends 
of friends method proposed by Huchra & Geller (1982) 
and the hierarchical clustering method (Materne, 1978; 
Tully, 1980, 1987). The latter method gave, on average, 
smaller groups than the former. The adopted final cata- 
log of groups was defined as that one which includes only 
groups (as well as galaxy members) common to the two 
catalogs. Accordingly, 3381 objects are field galaxies. Of 
the remaining grouped objects, 2703 galaxies are members 
of 485 systems with at least three members; most of them 
are groups which contain less than 10 members, and 13 
systems have more than 20 members (among which the 
clusters Virgo, Centaurus, Perseus). 
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A somewhat different approach to the construction of 
a large and deep sample of optical galaxies with good 
completeness in redshift was followed by Santiago et al. 
(1995). Their "Optical Redshift Survey" is a redshift 
survey of an optical galaxy sample which contains 8266 
galaxies with known redshift. This sample consists of two 
largely overlapping subsamples (which arc limited in mag- 
nitude and in diameter, respectively) drawn from the cat- 
alogs UGC, ESO, ESGC. The authors selected their own 
sample according to the raw (observed) magnitudes and 
diameters (rather than the corresponding corrected val- 
ues) and then quantified the effects of Galactic extinction 
(as well as random and systematic errors) on the optical 
density field. Redshift-distances were used in general, but 
the members of rich clusters were placed at each cluster's 
center. Probing the density field at large depth (up to 
cz — 8000 km s~^ ) requires much care in the deriva- 
tion of appropriate selection functions, which quantify the 
loss of galaxies due to magnitude or diameter limits or 
other observational selection effects. Taking into account 
possible non-uniformities in galaxy sampling, the authors 
attempted to match the original three catalogs from which 
their sample was drawn, by using separate selection func- 
tions (Santiago ct al., 1996). 

The reliability of the distances given in all above- 
mentioned galaxy samples (and hence of the resulting 
galaxy density field) is weakened by the fact that simple 
redshift-distances are in general adopted with no correc- 
tions (or no adequate corrections) for the peculiar motions 
which are due to inhomogeneous mass distribution and 
which prevent redshift to scale linearly with distance. 

In the present paper we propose to recover the true 
distances of a sample of nearby galaxies by treating ad- 
equately the effects of peculiar velocities of galaxies. As 
shown in a forthcoming paper (Marinoni et al., 1998b), 
in whic;h the resulting 3D galaxy density field is presented 
and local galaxy density parameters (corresponding to var- 
ious smoothing scales) are discussed, corrections of galaxy 
distances for peculiar motions appear to have a large im- 
pact on the values of the aforementioned parameters on 
small scales (i.e. on < 1 Mpc scale) and, hence, on the 
characterization of the environment (see §6). 

For our purpose, we take the above-mentioned complete 
sample extracted from LEDA (Garcia et al., 1993; Garcia, 
1993), for which all galaxy data as well as group assign- 
ments are available in the literature. Thus, we limit our- 
selves to the nearby universe (within cz = 5500 km s~^ ), 
where the peculiar velocity field is best known and no large 
corrections for galaxy sample incompletion are required. 
However, the galaxy sample considered is deep enough to 
cover interesting regions of prominent overdensities, such 
as the "Great Attractor" (GA) region and the Perseus- 
Pisces supercluster. 

The study of the peculiar velocity field resulting from 
a proper analysis of large redshift-distance samples (i.e., 
samples of galaxies having both redshift and redshift- 
independent distance estimates) has acquired a well- 
established role in the cosmological context (e.g., for test- 



ing the gravitational instability paradigm for the origin of 
large-scale structure, for deducing the relative distribu- 
tions of luminous and dark matter, and for constraining 
the value of the cosmological density parameter flo; see, 
e.g., the reviews by Dekel, 1994, and Strauss & Willick, 
1995). 

In general, peculiar velocity analyses proceed along two 
main lines of research. One can turn a distribution of 
galaxies or galaxy systems with known redshifts into a 
smoothed mass-density field (under the linear approxima- 
tion of the gravitational instability theory) and then pre- 
dict from it the peculiar velocity field , e.g., by means 
of the " dipole analysis" . Despite the complexity of these 
studies, the velocity vector of the Local Group (LG) with 
respect to the reference frame defined by the cosmic mi- 
crowave backgroimd (CMB) {vlg ^ 630 km s"-'^ towards 
I ~ 277° and b - 30° (Lubin & Villela, 1986; Kogut et al., 
1993)) has been approximately reconstructed in terms of 
the gravitational acceleration exerted by various popula- 
tions of cxtragalactic objects (e.g., the references cited in 
Kolokotronis et al., 1996). 

Alternatively, one can face the problem of deriving the 
true matter distribution from the observed peculiar veloc- 
ity field, e.g using few-parameter toy models or the non- 
parametric POTENT method developed by Bertschinger 
& Dekel (1989). 

In this paper, following the second kind of approach, we 
attempt to correct the Hubble flow for peculiar motions by 
means of two basic independent approaches to modeling 
the velocity field in the nearby universe: i) a modified clus- 
ter dipole model, which is Branchini & Plionis' (1996, here- 
after BP96) optical cluster dipole reconstruction scheme 
that we modify by including a local model of the Vir- 
gocentric infall in the LS region. We regard the cluster 
dipole model as a typical approach in which the peculiar 
velocity field is self-consistently derived through iteration 
techniques from the redshifts and positions of a sample of 
objects, ii) a multi-attractor model, which we take as an 
example of a parametric toy model which is fitted to cata- 
logs of peculiar velocities. In the multi-attractor model we 
assume that the velocity field is generated by a few promi- 
nent gravitational sources (Virgo cluster , Great Attractor, 
Perseus-Pisces Supercluster, Shapley concentration). 

We construct these models relying on the various data 
subsets contained in the Mark II and Mark III catalogs of 
peculiar velocities of galaxies. In particular, the Mark III 
catalog (Willick et al., 1997a), which is the most compre- 
hensive homogenized catalog of peculiar velocities avail- 
able today, is a merging of several data sets of spiral and 
elliptical galaxies. 

Through these peculiar velocity field models we provide 
homogeneous estimates of distances for our galaxy sam- 
ple. Moreover, we can evaluate to what extent differences 
in current views on the peculiar motions affect the recover- 
ing of galaxy distances and, hence, the 3D maps of galaxy 
distribution. 

The outline of our paper is as follows. In §2 we describe 
the properties of the Mark II and Mark HI catalogs which 
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are relevant to our velocity field models. In §3 we present 
our two basic approaches to modeling the peculiar velocity 
field and we discuss the results relative to different data 
sets. In §4 we address the inversion of the redshift-distance 
relation and, in particular, the ambiguity inherent to the 
triple-valued zones of this relation. We attempt to solve 
this problem using blue TF relations calibrated on suitably 
defined samples of galaxies having distances predicted by 
velocity field models. In §5, outlining some basic differ- 
ences among the velocity field models, we discuss the re- 
sulting sets of galaxy distances. In our conclusions (see §6) 
we mention some developments of peculiar velocity studies 
and focus on the effects of our galaxy distance corrections 
on the determination of the 3D galaxy density field. 

Throughout, the Hubble constant is 75 h km Mpc~^ 
In general, distances are expressed in velocity units 
(km ). 

2 The Mark II and Mark III peculiar velocity 

CATALOGS 

Our modeling of the peculiar velocity field is substan- 
tially based on several data subsets contained in the Mark 
III catalog of radial peculiar velocities (Willick et al., 1995, 
1996, 1997a). This is a homogenized database of redshift- 
independent distance estimates obtained using the revised 
Fabcr-Jacson (_D„ — a) relation for early-type galaxies and 
both the forward and inverse TF relations for spirals. In- 
put data come from different samples of galaxies, but the 
inferred galaxy distances arc the result of an accurate anal- 
ysis which carefully takes into account differences in the 
observational selection criteria, in the methods of mea- 
surement, in the TF calibration techniques, in statisti- 
cal bias effects. The uniformity of corrected data is en- 
sured through transformations onto a common system and 
through a mutually consistent TF calibration for the sam- 
ples of spirals. The Dn — a- inferred distances for early- 
type galaxies are rescaled in order to match the distance 
scale of spirals. 

This catalog is the current evolution of the previous 
Mark II catalog compiled by Burstein (1989), which is a 
merged set of 1184 galaxies grouped into 704 objects (sin- 
gle galaxies, groups, and clusters). Of these galaxies, the 
ellipticals (E) and lenticulars (SO) come from a combina- 
tions of the datasct used by LyndcU-BcU et al. (1988) in 
their analysis of the GA region and those by other au- 
thors (Lucey & Carter, 1988; Faber et al., 1989; Dressier 
& Faber, 1990) The spiral sample consists of the Aaron- 
son "good" and "fair" field spirals (see Faber & Burstein, 
1988, and Aaronson et al., 1989, for the definitions), the 
Aaronson cluster spirals, the spirals of dc Vaucouleurs & 
Peters (1984). In some cases we present results derived 
from the Mark II data alone in order to emphasize the 
differences between Mark II and Mark III data. 

With nearly 3400 galaxies, the Mark III catalog, which 
includes Mark II data with some improvements, is almost 
three times bigger than its predecessor and has a better 
space coverage. It contains new samples of spirals which 



cover sky regions poorly mapped in Mark II. Maps of the 
spatial distribution of the various subsamples, which probe 
different regions of the sky, are presented by Kollatt et al. 
(1996). In our analysis we use the following different sub- 
samples, avoiding the possibility that the same object be 
included in different samples: 

1. the E/SO galaxy sample (Lynden-Bell et al. 1988; 

Faber et al. 1989; Lucey & Carter, 1988; Dressier & 
Faber, 1990). It includes 544 galaxies grouped into 
249 objects which map the GA region. 

2. the HMCL cluster sample of spirals (Han & Mould, 
1992) with 36 objects. 

3. the WCF group sample of spirals (Willick, 1991; 
Courteau, 1992; Courteau et al., 1993; Courteau, 
1996) which contains 65 groups of galaxies covering 
the Perseus-Pisces region. 

4. the MAT data set of spirals (Mathewson, Ford & Bu- 
chorn, 1992) with 277 groups selected in the southern 
hemisphere. 

5. the A82 spiral sample (Aaronson et al., 1982) as re- 
vised by Tormen & Burstein (1995). Of the 359 
galaxies contained in the sample. 222 have been 
grouped in 67 groups (Willick et al., 1996) and 137 
are field galaxies. 

The analysis of different data subsets allows us to test 
the stability of our final results against the particular sub- 
sample used. The whole data set, in which we take to- 
gether spirals and ellipticals (which should trace the same 
velocity field) will help us to give more stringent con- 
straints on our models than smaller subsamples would. 

The accuracy of the peculiar velocities achieved in Mark 
III catalog is not only due to a careful homogenization of 
different data subsets and DIs for getting a uniform dis- 
tance scale; it is also due to the new approach with which 
the authors treat sample selection (or calibration) effects, 
which typically enter in the calibration of the DI being 
used in flux-limited samples, and the inferred-distance ef- 
fects, also called Malmquist bias, which typically enter 
when a calibrated DI is used to infer distances and hence 
peculiar velocities (e.g. WiUick, 1994). Selection effects 
are due to the fact that a magnitude limit in the selec- 
tion of the sample used for calibration at a fixed true dis- 
tance (e.g., in a cluster) tilts the direct TF regression line 
{Mb versus \ogVm) towards bright Mb at small values 
of logl^. The homogeneous part of the Malmquist bias 
arises from the space geometry in the sense that the in- 
ferred distance d underestimates the true distance r be- 
cause it is more likely to have been scattered by errors 
from r > d rather than from r < d, the volume being 
oc r^. The inhomogeneous part of the Malmquist bias 
arises from number density fluctuations and tends to en- 
hance systematically the inferred density perturbations. 

Clearly, removing the overall Malmquist effect (known 
as inhomogeneous Malmquist bias), arising from both vol- 



4 



ume effects and density variations, is of fundamental im- 
portance to yield unbiased peculiar velocities. In the Mark 
III catalog this was done using smoothed density fields ob- 
tained from the IRAS 1.2 Jy redshift survey (Fisher et al., 
1995), with the effect of peculiar velocities corrected for 
using linear theory (Yahil et al., 1991). The resulting field 
is expected to represent the general spiral density field on 
large scales. 

Thus, besides the uncertainties inherent to DI cali- 
bration procedures, the Mark III distances rely also on 
a model-dependent reconstruction of the general density 
field from redshift data. Moreover, redshift limits in a 
sample modify the nature of the Malmquist bias correc- 
tions applied. It would be worse to neglect any correction 
for the Malmquist bias; but we have to take into account 
the possibility that further errors and systematic effects 
can affect the data. Therefore, in our analysis we increase 
the tabulated errors on distances by a factor of 3, accord- 
ing to D. Burstcin's advices for the Mark II catalog. 

For all the single objects we have taken the Malmquist- 
corrected forward TF distances, since this is the best ap- 
proach to follow in our statistical analysis. Our approach 
belongs to the method I of the "method matrix" of pecu- 
liar velocity analysis illustrated in §5.2 of Strauss & Willick 
(1995). On the other hand, the forward TF distances of 
the spiral groups are fully corrected only for selection bias, 
because they are formed using redshift-space criteria and 
in this case it is selection bias rather than Malmquist bias 
which pertains. Residual Malmquist bias due to clustered 
structures still affects the data of spiral groups. Neverthe- 
less, the grouping procedure reduce distance uncertainties 
by a factor N'^/'^, if N is the number of grouped galaxies. 

3 Correcting Galaxy Distances through 
Velocity Field Models 

At the low redshifts discussed in this paper the Hubble 
law, cz = Hor,]s an excellent approximation in describing 
the Fricdmann expansion of the universe. However, galax- 
ies have motions above and beyond their Hubble reces- 
sion velocities. Therefore, by combining the cosmological 
component of redshift with the peculiar one, the Hubble 
relation is modified to: 

cz « H^r + [v(r) - v(0)]-f + + 

_ (v(0)-f)(v(r)-r) ^^^^ _ (v(r)-f)(v(0)-f) j^^^ 

where f is the unit vector towards the galaxy in ques- 
tion, v(r) is the peculiar velocity of a galaxy at position 
r, and v(0) is the non-comoving velocity of the observer. 

To first order approximation and in the Local Group 
rest frame, the distance of a galax;y (expressed in km s~^ ) 
can be written as: 

r = cz- [v(r,;,6) - v(0)]-f (2) 
where I and h are the galactic longitude and latitude. 



Thus, with a peculiar velocity field model which de- 
scribes point by point the function v(r, I, b), it is possible 
to correct the linear Hubble law in order to take into ac- 
count peculiar motions. We rely on two basic approaches 
to modeling the peculiar velocity field. They are described 
below. 

3.1 The Modified Cluster Dipole Model 
3.1.1 The Model 

We use the BP96 optical cluster dipole model as a basic 
predictive model of the peculiar velocity field. The BP96 
reconstruction scheme of the 3D positions and peculiar ve- 
locities of galaxy clusters is based on the observed distribu- 
tion of an optical cluster sample (with r < 25000 km s~^ ) 
extracted from the Abell/ACO sample (Abell, 1958; Abell, 
Corwin & Olowin, 1989). Taking into account observa- 
tional biases through an homogenization procedure and 
filling artificially the zone of avoidance with a simulated 
population of galaxies, BP96 obtained a statistically ho- 
mogeneous all-sky sample of clusters. 

Under some general and reasonable assumptions, e.g. 
that the peculiar velocities are caused by gravitational in- 
stability, that linear instability theory applies and that 
cluster density fluctuations are related to the mass den- 
sity fluctuations by a constant linear biasing factor 5c{v) = 
hc6{v), a simple relation between the pecuhar velocity of 
a galaxy and the surrounding mass density fleld can be 
derived (cf Peebles, 1980): 

where 3{r) = [p{r) — pb]/ Pb is the density fluctuation 
about the mean background density pb and (3 = ^^^^ ~ 
£^ 

Within the linear theory approximation, BP96 used an 
iterative reconstruction algorithm to remove redshift space 
distortions, recovering both the real space distribution of 
clusters and the fleld of peculiar velocities generated by 
this cluster distribution. 

The resulting density and velocity fields clearly depend 
on an assumed value of /3, which is not directly provided 
by the procedure itself. Nevetheless, this kind of analysis 
yields a resulting cluster 3D dipole pointing only ~ 10° 
away from the CMB apex; therefore, comparing the re- 
constructed cluster dipole amplitude to the LG peculiar 
velocity (as inferred by the CMB dipole) BP96 estimated 
13 ~0.21. 

The model gives some important information (such as 
bulk flow amplitudes) on large scales, but we are inter- 
ested in constructing an accurate picture of the velocity 
fleld and in recovering galaxy distances in the region which 
lies within about 5500 km s~^ . 

Because of the heavy smoothing procedure used, the 
cluster density is not known with a good local resolution. 
In other words, local contributions to the velocity fleld 
such as those originating from the gravity of the Virgo 
cluster, which is not even included in the Abell/ACO 
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Fig. 1. — The plots show the radial components (in the CMB frame) of the smoothed velocity fields in the supergalactic plane 
SGX, SGY for the cluster dipole model by Branchini & Plionis (1996) (left) and for the cluster dipole model modified with the 
inclusion of a Virgocentric infall (right). The field is smoothed with a 1500 km s~^ Gaussian window and is normalized to /3 = 0.21. 
The arrows and the boldface arrows distinguish between incoming and outcoming objects. The contours correspond to the same 
radial peculiar velocity; contour spacing is 100 km s~^ , with the heavy contour marking 0. Regions of different radial peculiar 
velocities are also indicated by shading. The Local Group is at the center, the Great Attractor is on the left, Perseus-Pisces is on 
the right, and Coma is at the top. Coordinates are expressed in km s~^ . 



cluster sample because of its large angular size, are not 
taken into account. 

We have address this problem by merging the cluster 
dipole model with a simple model of the Virgocentric infall. 
In our model we assume that the Virgo mass is spherically 
symmetric and is radially distributed with a King density 
profile 



p{r) = pb 



l + A[l 



-1.5 



(4) 



where x = r/rc, r is the distance of a galaxy from the Virgo 
center, rc is a smoothing length, and ^ is a normalization 
mass factor. Using eq. we obtain 



v(r) 



-AVI 



0.6„ 



In I 



1 



(5) 

where f is the unit vector pointing outward Virgo; the 
parameter plays the role of a smoothing length neces- 
sary for the linear theory to hold. 

We evaluate the influence of the Virgo cluster on the dy- 
namics of the LS by fitting the predictions of our modified 
cluster dipole model with the observed radial velocities of 



a sample of nearby galaxies having redshift-independent 
distances. 

In this fit we fix the coordinates of Virgo cluster at 
I = 284°, h = 74°, r = 1350 km s^^ (Virgo distance) as 
in Faber & Burstein (1988), whereas we treat C = AVt^-^ 
and Tc as free parameters of the model. Necessary ingredi- 
ents in the fitting procedure are the radial components of 
peculiar velocities as predicted by the model (?;™°'*), the 
corresponding observed recession velocities {czi) and the 
galaxy distances . The model parameters are determined 
by means of a minimization procedure; we minimize the 
quantity 



X 



(6) 



where cz™°'^ is given by eq. ^ and depends on a number 
of model parameters ak ■ The quantity is the quadratic 
sum of the uncertainties in the redshift measures, errors 
in the predicted redshifts induced by distance uncertain- 
ties plus a term which describes noise in the velocity field, 
under the assumption that these errors are independent. 
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Fig. 2. — The plots show the velocity field in the CMB frame for the cluster dipole model (left) and the modified cluster dipole 
model (right). The vectors shown are projections of the 3D velocity field in the supergalactic plane SGX, SGY . The contours 
correspond to the same velocity vector modulus; contour spacing is 100 km s~^ , with the heavy contour marking 100 km s~^ . 
Regions of different peculiar velocities are also indicated by shading. The Local Group is at the center, the Great Attractor is on 
the left, Perseus-Pisces is on the right, and Coma is at the top. Coordinates are expressed in km s~^ . 



The last term is a value of dispersion associated to non- 
linear and non-spherically symmetric motions of galaxies 
which lie around collapsed systems and for which our mod- 
els fail to predict the true radial velocity. For the first term 
we assume the value of 50 km s^^ , for the last we adopt 
the value of 200 km s^^ (e.g., Strauss & Willick 1995). 
We estimate the formal Icr errors of the model parameters 
from the covariance matrix calculated at the point which 
minimizes eq. |^ in parameter space. 

For this fit we choose the A82 spiral sample described in 
§2. These galaxies are distributed quite uniformily in space 
around the Virgo cluster and therefore they are suitable 
for studying the Virgo infall in detail. The forward TF 
relation used for estimating the distances of these galaxies 
has a rras scatter of cr = 0.47 mag (Willick et al., 1997a). 
Specifically, we consider the 158 galaxies of the A82 sam- 
ple which have distances smaller than 3000 km s^^ . This 
upper limit reflects the vanishing of Virgo gravitational 
influence at large distances. Moreover, the A82 sample is 
strongly incomplete at large distances, since it exhibits an 
abrupt reduction in the number of objects per unit redshift 
at czq = 3000 km s~^ . This redshift limit affects the re- 
liability of the inhomogeneous Malmquist bias corrections 



for the 59 objects lying beyond that limit. 

3.1.2 Model- fitting results 

The resulting best-fitting parameters C — AVp-^ and Tc 
(km s~^ ) are tabulated in Table 1 together with the value 
of per degree of freedom, which measures the good- 
ness of the fit, and the rms dispersion of model velocity 
residuals cr{cz) (in km s^"'^ ). 

TABLE 1 

The best-fitting parameters for the 
modified cluster dipole model. 

C rc y?ldof <r(cz) 

(km ) (km s"-*^ ) 

3.4 ±1.3 810 ± 110 097 404 



Fig. 1 shows the smoothed radial velocity fields (in 
the supergalactic plane SGX, SGY) of the BP96 cluster 
dipole model and our modified cluster dipole model. Fig. 
2 shows the smoothed velocity fields in the supergalactic 
plane SGX, SGY for the same models. 

Although the BP96 cluster dipole model suffers from 
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uncertainties in the relation between fluctuations in mat- 
ter and light distributions (e.g the linear biasing assump- 
tion), from sample selection effects, and from the method 
used to fill artificially the Zone of Avoidance, the model 
fits fairly well the A82 data in the LS region. The 
model implies a scatter of a{cz) = 485 km s^^ (with 
X^/dof = 1.44) which can be compared with the value 
of cr(cz) — 554 km s~^ (with x^/dof = 2.14) which one 
would obtain by fitting an unperturbed Hubble flow (in 
the CMB frame) to the same data. 

The fit becomes even better if we add the perturbing 
presence of Virgo cluster (see the low values of cr(cz) and 
X^/dof listed in Table 1). From our best fit we infer that 
the LG has a Virgo infall velocity equal to v = 440 ± 180. 
This value is larger than those quoted recently by several 
authors for whom the Virgo infall is not a major source of 
the velocity field in the LS, in disagreement with earlier 
contentions (see Han & Mould, 1990 for a summary table). 
In recent years, a large Virgo infall amplitude was also sug- 
gested by Tonry et al. (1992), who, applying the surface 
brightness fluctuation technique as DI for local ellipticals, 
found the value of v — 340 ± 80 km s^^ . On the other 
hand, from the analysis of the magnitudes of bright clus- 
ter galaxies in a way which is free of assumptions about 
motions in the LS, Gudchus (1995) found no signiflcant 
evidence of a Virgo infall of the LG. 

The too great an amplitude of LG infall velocity is not 
due to differences in the choice of observed Virgo redshift. 
It is instead likely due to the fact that the BP96 cluster 
dipole model generates on the LG a velocity vector (with 
projections along the supergalactic axes vsgx — — 225±61 
km s-i , vsGY = 131 ± 40 km s'^ , vsgz = -349 ± 81 
km s~^ ) which has a lower component in the general 
direction of GA {vlg = 285 km s^^ ) than that pre- 
dicted by other studies {v ^ 500 km s^^ ; see, e.g., Faber 
& Burstein, 1988). So, in attempting to reproduce the 
local observed flow pattern, our fitting parameters have 
to compensate the underestimated GA infall of LG. This 
interpretation is consistent with the fact that the CMB 
anisotropy dipole has a component directed towards the 
Virgo cluster of u = 418 ± 25 km s~^ (see, e. g., the re- 
view by Davis & Peebles, 1983) and is supported by the 
fact that, by adding the Virgo mass contribution, we find 
a velocity component of the LG motion towards GA of 
Vlg — 558 ± 170 km s^^ , in good agreement with the 
value of Vlg — 535 km s^^ obtained by Faber & Burstein 
(1988) with an independent model. 

3.2 The Multi-Attractor model 

3.2.1 The Model 

The GA model of the Seven Samurai, proposed by 
Lynden-Bell et al. (1988) and revised by Faber & Burstein 
(1988), is the current most popular phenomenological 
model of the large scale peculiar motions in the nearby 
universe. As claimed by several authors in the last years 
(e.g., Faber & Burstein, 1988; Han & Mould, 1990; Shi- 
masaku & Okamura, 1992) and substantially confirmed by 



a look at the mass density fields reconstructed through the 
application of the POTENT method to Mark II and Mark 
HI data (e.g., Dekel et al. (1993) and Sigad et al. (1998), 
respectively), most of the local velocity field can be in- 
terpreted as being generated by few gravitational sources 
characterized by spherical symmetry. It is then sensible to 
describe the velocity field in term of a multi-attractor toy 
model dominated by spherically symmetric Virgocentric 
and GA-infalls. In this model we consider also the possible 
effect of other gravitational sources, such as the Perseus- 
Pisces Supercluster and the Shapley concentration. 

This kind of model is an oversimplified version of facts, 
because it imposes a simple, spherical geometry on the 
adopted gravitational sources of a complex velocity field, 
which actually arises from a continuous field of asymmet- 
ric density fluctuations. Nevertheless, this model, which 
as yet has never been applied to Mark HI data, is the 
most simple and statistically significant tool to analytically 
correct the Hubble law. As discussed below, the multi- 
attractor model yields predictions in satisfactory agree- 
ment with data and it is adequate for our purposes. 

Instead of considering the gravitational instability 
paradigm, we use the infall model discussed by Regos & 
Geller (1989). In this model, based on the Friedmann solu- 
tion, the motion of galaxies at a particular radius from the 
cluster center is approximated as the motion of a spherical 
mass shell, which follows the same gravitational equation 
of motion as the expansion factor of the universe. The 
shell can be treated as a Friedmann universe on its own; it 
is possible to define a formal Hubble constant and a formal 
cosmological density parameter for the shell. This infall 
model yields an implicit exact dependence, and an explicit 
approximate analytic expansion, for the peculiar velocity 
as a separable function of the present density contrast (^o) i 
averaged inside a radius r, and the present cosmological 
parameter ilo- In the range 0.1 < JIq < 1 this espansion 
can be written (Regos & Geller, 1989) as: 



O.06317;]'^«((5o 



0.027r!O-^i((5o)^ - 0.015n°-^^{So)* 



o.omr'('5o)^ 



(7) 



where the distance r is expressed in km . 

This series expansion provides convergence to the exact 
solution for (Sq) < 1 — 2. However because the linear term 
dominates in the peculiar velocity, the Qq dependence of 
the exact solution is close to that of the linear approxima- 
tion. 

We approximate the total peculiar velocity as the sum 
of the components generated by the individual attractors. 
Assuming a King profile (eq. U) for each attractor (so that 
infall velocities converge at the center of attracting mass) , 
we build up a parameter toy model which is fitted to the 
various galaxy subsamples which are part of the Mark III 
data set. 
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Eq. is strictly applicable only in the case of a single 
attractor, as peculiar velocities can be added under the 
hypothesis, valid in linear theory (and even beyond, see 
Susperregi & Buchert, 1997) that they are proportional to 
peculiar acceleration. However, eq. ^ reduces to linear the- 
ory almost anywhere, except in the proximity of an attrac- 
tor, where the contributions from all the other attractors 
are weak. Therefore, the non-linear formula is considered 
as a suitable approximation for the velocity field generated 
by the attractors. As a consistency check, we have verified 
that the use of pure linear theory does not change signifi- 
cantly the obtained values of the free parameters, with the 
exception of flo, which can not be disentangled from the 
normalization mass factor A of eq. ^ if no non-linearity is 
taken into account. 

Specifically, we assume that the perturbing masses are: 
1) the Virgo cluster (V), which gives rise to distortions in 
the Hubble flow of LG surrounding regions (e.g., de Vau- 
couleurs & BoUingher, 1979; Tonry & Davis, 1981); 2) the 
Great Attractor (GA), which seems the best candidate to 
explain large scale galaxy motions inside a volume of ra- 
dius 5000-10000 km s'^ (Burstein et al., 1986; Dressier et 
al., 1987); 3) Perseus-Pisces (PP), a large and irregular su- 
percluster, whose massive filamentary structure seems to 
be concentrated at a redshift of about 5000 km s~^ in the 
antipodal direction of GA (Willick 1991; Han & Mould 
1992; Courteau et al., 1993); 4) the Shapley concentra- 
tion (SH), a very rich concentration of ~25 galaxy clus- 
ters whose mean complex lies at 14000 km and is lo- 
cated behind the GA region (Scaramella et al., 1989). SH 
stands out as the richest optical supercluster of the entire 
sky within z < 0.1 (Zucca et al., 1993). 

In our final analysis we do not consider the possible grav- 
itational influence of the Coma Supercluster, because we 
have checked that its inclusion does not improve the good- 
ness of the fit (see, e.g., Shaya, TuUy & Pierce, 1992, who 
reached the same conclusion in their less sensitive analysis 
of a galaxy sample limited to cz = 3000 km s~^ ). 

In the fit we do not consider the few very nearby galaxies 
having cz < 700 km and we do not correct their tabu- 
lated distances. These galaxies are believed to share with 
the LG a bulk motion (the so-called Local Anomaly) per- 
pendicular to the supergalactic plane (Faber & Burstein, 
1988). The Local Anomaly is a perturbation term intro- 
duced in the standard picture of the peculiar velocity field 
of the very nearby universe in order to explain the dis- 
crepancy between the observed and the predicted radial 
velocities of the LG with respect to the CMB. It is in gen- 
eral interpreted as a due to the " negative gravity" of the 
"Local Void" (TuUy & Fisher, 1987) which is located above 
the supergalactic plane. 

After we have specified for each galaxy the redshift czi 
and the distance , the expression (see eq. ||) becomes 
a known function of several parameters. The free param- 
eters which describe mass distribution and geometry are, 
for each attractor, Ai, Vd- Furthermore, in the case of GA 
and SH, we have decided to leave free the Galactic coordi- 
nates l,h of their centers and also the distance of GA from 



us. Finally, we try to leave flo as a further free parameter. 
We take, instead, from the literature the values for the po- 
sitions of PP {I = 120°, b = -30°, r = 5000 km s"^ ) and 
\ {1 = 284°, 6 = 74°, r = 1350 km s^^ ), and the distance 
of SH (r = 14000 km s^^ ). 

Thus, we are left with 13 degrees of freedom. Solu- 
tions may be unstable or physically not meaningful in the 
presence of so many parameters. However our data set 
is very large and quite uniform in spatial distribution so 
much to prevent statistical artifacts. Moreover, the im- 
proved accuracy in estimated distances ensures confidence 
in the possibility of constraining different models through 
the minimization procedure. Lastly, if we assume that the 
significance of each parameter is related to the change in 
likelihood of the best fit that results from its addition to 
the model, our trials with models having less parameters 
confirm the validity of the best fits achieved with many 
degrees of freedom. Nevertheless, the minimization 
scheme may introduce biased results due to the incom- 
pleteness of the dataset. For example, the redshift limit of 
Mark HI catalog, which does not map distant regions of 
the universe, could constrain in a wrong way the parame- 
ters of the Shapley concentration. 

The significance is ill-defined because velocity errors 
are coupled and because it depends considerably on the 
size of the distance errors and velocity spread due to noise 
in the velocity field. Encouragingly, we have checked that 
different values of the velocity spread (taken as constant 
in space) have little influence on the attractor parameters. 
In any case, we do not use the statistics to assess the 
validity of a model. We are interested in the relative de- 
crease of the x^ values, because the purpose of our analysis 
is to obtain a redshift-distance relation having less scatter 
and systematic bias than that of the simple unperturbed 
Hubble flow. 

3.2.2 Model- fitting results 

Using various subsamples we estimate the attractor pa- 
rameters separately in order to examine their reliability 
and stability. 

In Tables 2 and 3 we present the data number N , the 
values of x^ jdoj and rms dispersions of model velocity 
residuals (t(cz) (in km ) relative to Mark II and Mark 
HI data for Hubble flows in the CMB frame and in the LG 
frame, respectively. 

We report results for the whole Mark II and for a Mark 
II subsample (Mark II*) which docs not include the spirals 
by de Vaucouleurs & Peters (1984) and the field ellipticals. 
The omission of the former subsample is motivated by the 
fact that it docs not appear in the subsequent compilation 
(Mark HI). The omission of the latter is justified by the 
fact that many ellipticals reside in systems of fairly high 
velocity dispersion and hence, can not be retained as good 
tracers of the velocity field. In these tables we give also 
the results relative to the various Mark III subsamples (de- 
noted as described in §2), to the Mark HI spirals (Mark 
HI*) and the entire Mark III. 

In Table 4 we give the minimization results (together 
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Fig. 3. — The plots show the radial components (in the CMB frame) of the velocity fields in the supergalactic plane SGX, SGY 
for the multi-attractor model fitted to the Mark II subsamplc (left) and the whole Mark III (right). The arrows and the boldface 
arrows distinguish between incoming and outcoming objects. The contours correspond to the same radial peculiar velocity; contour 
spacing is 100 km s~^ , with the heavy contour marking 0. Regions of different radial peculiar velocities are also indicated by 
shading. The Local Group is at the center, the Great Attractor is on the left, and Perseus-Pisces is on bottom right. Goordinates 
are expressed in km s~^ . 



Fig. 4. — The plots show the velocity field in the CMB frame for the multi-attractor model fitted to the Mark II subsample (left) 
and the whole Mark III (right). The vectors shown are projections of the 3D velocity field in the supergalactic plane SGX, SGY . 
The contours correspond to the same velocity vector modulus; contour spacing is 100 km s~^ , with the heavy contours marking 
200 km s~^ and 300 km s~^ for the Mark II subsample and the Mark III, respectively. Regions of different peculiar velocities are 
also indicated by shading. The Local Group is at the center and the Great Attractor is on the left. Goordinates are expressed in 
km s-i . 10 



TABLE 2 

The Hubble flow model in the CMB frame. 



Mark II* Mark II A82 MAT HMCL WCF E/SO Mark III* Mark III 



xVdof 
cr{cz) 



413 


662 


192 


277 


36 


65 


250 


570 


820 


1.90 


1.70 


1.43 


1.16 


1.97 


1.30 


3.00 


1.22 


1.97 


698 


835 


1020 


634 


696 


851 


977 


812 


873 



TABLE 3 



The Hubble flow model in the Local Group frame. 





Mark 11* 


Mark II 


A82 


MAT 


HMCL 


WCF 


E/SO 


Mark III* 


Mark III 


M 


413 


662 


192 


277 


36 


65 


250 


570 


820 


xVdof 


2.10 


1.71 


1.22 


1.49 


1.96 


1.34 


3.24 


1.41 


1.97 


cr{cz) 


682 


790 


943 


661 


627 


1032 


955 


804 


858 



TABLE 4 



The BEST-FITTING PARAMETERS FOR THE MULTI-ATTRACTOR MODEL. 





Mark II* 


Mark II 


A82 


MAT 


HMCL 


WCF 


E/SO 


Mark III* 


Mark III 


M 


413 


662 


192 


277 


36 


65 


250 


570 


820 


xVdof 0.94 


0.99 


0.69 


1.08 


0.97 


0.78 


2.60 


0.98 


1.48 


aicz) 


543 


686 


757 


618 


415 


739 


903 


713 


787 


Ay 


3.1±0.6 


2.9±0.8 


4.5±0.8 











2.0±1.2 


4.0±0.95 


3.1±1.6 


Aga 


2.9±0.5 


3.0±0.7 


3.7±0.9 





3.2±1.0 


3.9±1.4 


3.5±0.7 


1.5±0.5 


2.4±1.0 


App 


1.0±0.7 




3.2±1.7 





4.0±1.2 


2.8±1.0 





3.4±0.7 


2.6±1.1 


Ash 


1.0±0.9 







5.0±2.0 


4.2±1.2 





2.0±0.9 


4.6±1.4 


3.3±1.4 


rc,v 


530±30 


500±30 


500±28 











490±44 


480±73 


498±30 


1"c,GA 


2080±45 


1740±70 


1950±82 





1830±70 1930±110 


2500±68 


1820±112 


1806±70 


rc,PP 


1480±150 




1440±110 





900±50 2200±110 





1410±60 


1415±61 


rc,SH 


2450±170 







2500±130 


3000±50 





2725±80 


2930±50 


3004±35 


Iga 


305° ±4° 


309°±4° 


294° ±6.5° 




314° ±3° 303° ±3° 


309° ±3° 


300° ±18° 


309° ±5° 


bGA 


15°±3° 


18° ±3° 


17° ±4° 




0°±6° 


5° ±3° 


10° ±3° 


25°±19° 


18°±3° 


dcA 


4170±40 


4170±50 4150±30 




4200±90 4200±40 


4400±110 


4200±105 


4200±30 




312°±9° 






319° ±3° 


319° ±6° 




308° ±4° 


308°±10° 


308° ±5° 




27° ±5° 






8°±5° 


3° ±4° 




2° ±5° 


2°±13° 


4°±5° 


Ho 


0.8±0.2 


1.0±0.2 


0.6±0.3 


0.5±0.3 


0.5±0.2 


0.8±0.3 


0.2±0.2 


0.6±0.2 


0.4±0.3 



with 1(7 errors) relative to Mark II and Mark III for our 
multi-attractor model; the values of the GA distance da a, 
attractor radii rc, and model velocity residuals a{cz) are 
given in km s~^ . Wc have checked that the inclusion of PP 
and SH did not improve significantly the fit to the whole 
Mark II, so that only two attractors (Virgo and GA) are 
adopted in this case. 

For some parameters, in Fig. 5 we show the ellipses of 
confidence (at the 68% and 90% significance levels) ob- 
tained from the covariancc matrix, according to the model 
fitted to the whole Mark III. Fig. 3 shows the radial ve- 
locity fields (in the supergalactic plane SGX, SGY) of the 



multi-attractor models relative to the Mark II subsample 
and the whole Mark III. Fig. 4 show the respective velocity 
fields. 

From an inspection of Tables 2,3, and 4 we can draw the 
following comments. Velocity field models improve the fit 
with respect to the unperturbed Hubble flow in a decisive 
way, and a multi-attractor model is the best way to trace 
out the underlying observed velocity field. Only for the 
MAT sample, a bulk motion generated by SH seems to be 
best solution. On the other hand, the fact that the imi- 
form expansion model in the LG frame generally fits the 
data much better then that in the CMB frame strengthens 
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the evidence for a large streaming motion occuring along 
the PP-GA baseline and flowing towards SH. 

We stress that for most Mark III samples the best fit 
is obtained if SH is included in the velocity field models. 
Only in A82 and WCF subsamples the long range action 
of this supercluster is undetected, probably because the 
former sample maps only the LS dynamics and the latter 
mostly covers the region surrounding PP, which introduces 
incompleteness biases in the statistics. 
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Fig. 5. — For some parameters, we show the regions of confi- 
dence (at the 68% and 90% significance levels) according to the 
model fitted to the whole Mark III. 

In every subsample the application of velocity correc- 
tions reduces considerably both scatter and x^- 

The reduction in scatter is more evident in the A82, 
HMCL, WCF subsamples which have more accurate dis- 
tance estimates (see the relative cr- values in Tables 2, 3, 
and 4). Figs. 6, 7, and 8 show the recession velocity 
(in the CMB frame) versus the distance before and af- 
ter correction for peculiar motions as predicted by models 
fitted to the three aforementioned subsamples. These fig- 
ures clearly illustrate how closer to the Hubble line our 
corrected data lie than original data do. 

An independent way to judge the validity of a model is 
to look at the velocity deviations from model predictions 
as a function of the distance. Fig. 9 shows that, for the at- 
tractor models of the Mark II subsample and whole Mark 
HI, the residuals {observed minus calculated velocity) ex- 
hibit no overall systematic variation with distance, as is 
expected for a good model. 

The position of the GA obtained from various Mark HI 
subsamples is almost similar and stable. The whole Mark 
II and Mark III yield I = 309° b = 18° for the GA as in 
Fabcr & Burstcin (1988). Our best fit of the whole Mark 
HI locates the mass center of SH at I — 308° 6 = 2°, 
which somewhat differs from the position of the optical 



center {I = 319°, b = 27°) as reported by Scaramella et 
al. (1989). Interestingly, all Mark III subsamples tend 
to locate SH at low ^-values with respect to the optical 
position; this suggests that there is more matter behind 
the Galactic plane than previously thought. Ongoing red- 
shift surveys in the SH region (e.g., Proust, Quintana & 
Slezak, 1998) will be able to investigate on the possible SH 
extension towards the zone of avoidance. 

The relative values of the radii Vc give some insight into 
the distribution of matter inside the attractors. The large 
fc-values of GA and SH indicate that they are not peaks 
of highly clumped matter, but rather large overdensity re- 
gions. In particular, the GA mass is substantially dis- 
tributed across a distance of ~ 2000 km s'^ (Faber & 
Burstein (1988) quoted the value of 1500 km s~^ for a 
different model); the PP core radius is ~ 1400 km s^^ ; 
the SH core radius is ~ 3000 km , whereas the op- 
tical radius, defined as the radius of a sphere centered in 
SH and containing ^25 (or ^15) rich clusters, is ^ 5000 
km s-i (or -1500 km s^^ ) (Vettolani ct al., 1990; BardcUi 
ct al., 1994). 

The normalizing mass parameters Ai, which arc left as 
free parameters, bear no immediate interpretation. How- 
ever, these values are used to evaluate interesting physical 
quantities, such as the masses of the attractors or the infall 
velocities that they generate. 

To estimate the excess mass in the attractor regions, 
we evaluate the volume integral of the mass distribution 
we adopted (the King profile of eq. (4)). This is a crude 
approximation because of the possible strong departure of 
the mass distribution from spherical symmetry; we sup- 
pose that anisotropies in the density field cancel out when 
averages over a large volume of space. For a sphere cen- 
tered in GA and having LG on its border, the whole Mark 
II yields a mass of (5.1 ± 2.5) • 10^^ • ft-^Mg, whereas 
Lyndell-Bell et al. (1988), using only the elliptical sample, 
reported the value of 5.4- IO^^Mq, and Shaya et al. (1992) 
reported 1.5 • (iEj ^^ ■ lO^^^Mg. The whole Mark HI gives a 
smaller mass of (1.8 ± 1.6) • 10^^ h''^ Mq. 

If we consider the attractors as spheres of radius rc, 
we obtain the excess masses reported in Table 5 for our 
multi-attractor models fitted to some data subsets (Mark 
II subset, Mark III spirals, whole Mark III). Table 5 lists 
also the density contrast averaged over a sphere centered 
on the; attractor and having the LG on its border. 

According to Mark HI and Mark II, Virgo is con- 
firmed to be a poor cluster; it has an excess mass of 
1.3 • 10"ft-^MQ (see Table 5) and a total mass of 2.0 • 
lO^/i-^M© according to Mark HI data (2.2 • lO^/i-^M© 
and 3.2 • 10^^h~^MQ according to Mark HI spiral data, re- 
spectively). These values are substantial agreement with 
the estimates of masses based on optical and X-ray data. 
Prom optical data Girardi et al. (1998) estimated a virial 
mass of M = (2.7lo4) lO^^/i^^Afo within a virialization 
radius of 1.69h-^ Mpc. From ROSAT PSPC X-ray ob- 
servations Nulsen & Bohringer (1995) estimated a cluster 
mass per unit length of 1.24 • 10^^ Mq kpc~^, which is 
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Fig. 6. — Those plots show the recession velocity cz (in the CMB frame) as a function of the distance r (in km/s) before (left) 
and after (right) corrections for peculiar motions, as predicted by the multi-attractor models fitted to the A82 galaxy sample. The 
scatter decreases from a = 1020 (left) km to a = 757 km (right). 
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Fig. 7. — The same as in Fig. 6 for the multi-attractor model fitted to the WCF sample. The scatter decreases from a = 851 

km s"^ (left) to cr = 739 km s"^ (right). 
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Fig. 8. — The same as in Fig. 6 for the multi-attractor model fitted to the HMCL sample. The scatter decreases from cr = 696 
km s"^ (left) to = 415 km s"^ (right). 



consistent with our estimates, if we use an isothermal 
sphere extrapolation (M(r) oc r). 

Compared to Mark II, Mark III data stress the role 
of PP and SH (sec also the ^j-vahies listed in Table 4), 
whilst they give a less prominent GA with a mass which 
is roughly comparable to the PP mass and is ^ 5 times 
smaller than the SH mass. This is in good agreement with 
Scaramella et al.'s (1989) finding that in optical maps the 
SH region contains 5 times as many rich clusters as the 
GA region does. The SH excess mass derived from Mark 
HI (3 - 5.7 • IO^^H-^Mq, see Table 5) (as well as the cor- 
responding SH total mass of 4.7 - 8.1 • h''^ ■ IO^^Mq) is 
roughly consistent with values intermediate between up- 
per and lower limits based on optical and X-ray data. 
From the velocity dispersion of the SH clusters Raychaud- 
haury et al. (1991) and Quintana et al. (1995) obtained 



upper limits to the total mass of 1.3 • lO^^h ^ Mq and 
7- lQ^^h~^ Mq, respectively. 

Assigning to the 40 SH clusters the virial mass of 
3 • W^^hr^ Mq, which is the average of the virial masses 
deduced for 10 SH clusters, Quintana et al. (1995) esti- 
mated a lower limit for the SH mass of 1.2 • 10^^h~^ Mq. 
From the X-ray luminosities of the 12 brightest clusters, 
scaling to the Coma cluster (with Lx oc M° **), Raychaud- 
hury et al. (1991) obtained a lower limit of 9- 10^^ H'^Mq. 
From the analysis of the ROSAT PSPC and Einstein Ob- 
servatory IPC X-ray observations Ettori, Fabian & White 
(1997) derived a lower limit for the dynamical SH mass 
of 6 • lO^'^h^^ Mq; their virial mass estimate is ^5 times 
larger than the last value. From the SH influence on local 
flows Shaya et al. (1992) reported a large dynamical mass 
of 2.6 • fig-^ 10^'^ Mq. The same authors obtained a PP 



TABLE 5 

The excess masses and density contrasts for the attractors. 

Mark III Mark Ill(spirals) Mark Il(subset) 

Attractors M(< Vc) (So)lg M(< rc) (So)lg M(< rc) {So)lg 

(Mq/i-I) (Mg/l-l) 

Virgo (1.3± 1.3) • 10^** 0.36 ± 0.22 (2.2 ± 1.5) ■ lO" 0.44 ±0.23 (2.8 ± 1.5) ■ lO" 0.35 ± 0.30 

GA (4.7 ± 4.4) • 10^" 0.38 ±0.18 (4.4 ± 3.0) ■ lOi'"' 0.24 ±0.11 (1.6 ± 0.6) • lO"' 0.44 ±0.10 

PP (2.4 ± 2.2) • IQi^^ 0.19 ±0.09 (4.7 ± 2.7) • lO^^^ 0.23 ±0.06 (1.5 ± 2.1) • 10^^ 0.07 ± 0.07 

Shapley (3.0 ± 2.7) ■ IQi^ 0.12 ± 0.05 (5.7 ± 3.7) • IQi^ 0.16 ±0.06 (8.0 ± 9.8) • lO^^ 0.02 ± 0.02 
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Fig. 9. — We show the residual {observed minus calculated) ve- 
locities versus the distance r (in km s~^ ), for the velocity mod- 
els fitted to the whole Mark III (a) and the Mark II subsample 
(6). 

dynamical mass of 5 • Og* 10^^ Mq, which better agrees 

with our estimates than their SH mass does. The best- 
fitting values of Q,o tabulated in Table 4 are used for the 
estimates given above. 

We have also computed the infall velocity generated by 
Virgo alone at the LG position. We find v = 170 ± 64 
km s~^ , V = 145 ±57 km s~^ , and t; = 96±45 km for 
the A82 sample, the Mark III spiral sample, and the whole 
Mark III, respectively. In particular the last low value indi- 
cates that Virgo infall is not a major source of the velocity 
field in the LS, in agreement with several recent results 
and in disagreement with the results of the modified clus- 
ter dipole model (see end of §3.1.2). 

Remarkably, the influence of the GA at the LG position 
is lower than expected from Mark II and from previous 
results. 

TABLE 6 

The motion of the Local Group in the 
cmb frame. 

Sample 1 h 9 V V]^ 

Mark Il(subset) 305° 25° 25° 705 567 

Mark II 307° 29° 26° 666 563 

A82 290° 30° 11° 479 615 

MAT 310° 7° 48° 412 491 

HMCL 318° 0° 49° 810 410 

WCF 310° 0° 43° 225 456 

E/SO 309° 12° 34° 556 516 

Mark Ill(spirals) 309° 7° 37° 454 496 

Mark III 309° 12° 36° 488 511 



The pull generated by GA at the LG position turns out 
to be i; = 314 ± 200 km s~^ according to Mark III data 
(it is t; = 620 ± 260 according to Mark II data). There- 
fore, nearly one half of the LG motion is due to the SH 
action. At variance with Mark II, Mark III data suggest 
that the SH, projected behind the GA, yields a pull on the 
LG comparable to that exerted by GA. 

There is a way to test the relative prominence of the 
GA without resorting to any specific attractor modeling. 
It suffices to assume that the infall on the attractor be ra- 
dial and that we know fairly well its angular position in the 
sky. Let us take a sphere centered at half the distance be- 
tween LG and GA, with LG on its border. Now if the infall 
is radial and GA is the only responsible for the observed 
flow patterns, then, for a simple geometrical reason, the 
galaxies lying inside the sphere should have positive pe- 
culiar velocities, whilst the galaxies lying outside should 
have negative radial components (see Fig. 10); only at the 
surface of the sphere the radial components would become 
null. 










Great 
Attractor 


dc /' 








Local 
Group 



Fig. 10. — The plot illustrates the following concept: for a 
sphere centered at half the distance between the Local Group 
and the Great Attractor, under the hypothesis of a radial in- 
fall on the Great Attractor only, a galaxy which is at distance 
dc from the center of the sphere and lies inside (outside) the 
sphere should have a positive (negative) radial peculiar velocity 
Vp (as viewed from an observer placed in the Local Group). 

In Fig. 11 we compare the results obtained with Mark II 
and Mark III data. Mark II and A82 data (see Fig. lid, f, 
respectively) favour a massive GA dominating large scale 
motions, whilst evidence for this is considerably weakened 
in other Mark HI samples, which do not show a clear radial 
infall pattern towards GA. This can be only partially due 
to the perturbing presence of PP, because if we exclude 
the galaxies located in the PP region from the plot (see 



15 



2000 
1000 


-1000 
-2000 
2000 
1000 
^ 
^ -1000 
-2000 
2000 
1000 


-1000 
-2000 



-yy-.c-y/A 



I , ,h:r'T''r'r 



u-^':y.y::-: 



■:6^yy 




.'i'.'.r.'.'i.'.'.l1 I L 



^•.f , % ^ H 




2000 4000 6000 8000 



2000 4000 6000 8000 10000 



d, (km/s) 



Fig. 11. — The observed radial peculiar velocities Vp (in km s~^ ) as a function of the distance dc (in km s^^ ) from the middle 
point of a line joining the Great Attractor with the Local Group. If the Great Attractor were the only source of radial infall, then 
for dc > 2100 km s^^ peculiar velocities should be negative, while for dc < 2100 km s^^ they should be positive (see Fig. 10). We 
consider various samples of galaxies: the whole Mark III data (o), the Mark III spirals (fe) (almost all Mark III ellipticals lie in the 
GA region), the MAT 82 spirals (c) , the A82 spirals (d). the Mark III sample without the objects located in the Perseus-Pisces 
region (i.e., those placed inside a sphere centered on this supercluster and having a radius of 3000 km .) (e), the whole Mark 
II data (/). 



Fig. lie), there are still many galaxies with positive 
peculiar velocities outside the sphere. 

This reflects the lack of a back-side infall towards GA or, 
in other words, the long range effect of SH which prevents 
galaxies located in the back of GA to acquire negative pe- 



culiar velocity components. 

If the GA were the dominant system we would expect 
a GA back-infall to be comparable in amplitude to the 
GA forward-infall. On the contrary, if SH exerts a signif- 
icant gravitational influence, we expect an asymmetry in 
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the amplitudes of the back- and forward-infalls, i.e. pos- 
itive components outside the sphere. In conclusion, local 
samples such as the A82 sample and Mark II, which does 
not include many galaxies in the back of GA, can not fully 
reveal the importance of SH in determining peculiar mo- 
tions. 

In Fig. 12 we show the mean overdensity profile in- 
side a spherical volume centered on each attractor, for our 
model of the whole Mark III. The PP and GA overden- 
sities at r ~ 1200 km s""'^ are in reasonable agreement 
(being just a bit higher) with the density peaks at 1200 
km s""'^ Gaussian smoothing, {5o ~1.2 and ^1.4, respec- 
tively) found by Sigad et al. (1998) who used the PO- 
TENT procedure to reconstruct the smoothed mass den- 
sity field from Mark III data (for i}o=l), taking into ac- 
count mild non-linear effects. These mass density peaks 
are not much different from the reconstructed real-space 
density peaks (at 1200 km s^^ Gaussian smoothing) of 
IRAS galaxies (Sq ^0.8; sec Sigad ct al., 1998) and opti- 
cal galaxies (Sq ^1.8 and 0.8, respectively; see Hudson ct 
al., 1995), which is in line with the widespread contention 
that mass and galaxies are related via an approximate lin- 
ear biasing relation with a biasing factor be of order unity. 

In Table 6 we report the direction {I, b) and the am- 
plitude V of the CMB anisotropy dipole generated by 
the assumed distribution of attracting masses, for multi- 
attractor models relative to various subsamples. We also 
give the angular separation 6 between the observed and re- 
constructed dipole together with the projection Vj| of the 
observed CMB dipole in the {I, b) direction. In all cases the 
reconstructed dipole shows a satisfactory agreement both 
in direction and in amplitude with the observed dipole. 
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Fig. 12. — The plots show the density contrast, averaged over 
spheres centered on each attractor, as a function of the distance 
r (in km s~^ ) from the attractor, for the velocity model fitted 
to the whole Mark III data. 



This can be taken as an independent evidence of the 
reliability of the multi-attractor toy model we used. 

The remaining differences between the observed and re- 
constructed dipole can be justified in several ways. For ex- 
ample, inhomogeneities beyond 15000 km s^^ could con- 
tribuite significantly to the LG motion (e.g., Branchini, 
Plionis & Sciama, 1996; Plionis & Kolokotronis, 1998) or, 
according to an alternative explanation, the motions of 
the very nearby galaxies are also subjected to the Local 
Anomaly (Burstein, 1990). If we preferred the latter ex- 
planation, we would reconstruct a Local Anomaly velocity 
vector with an amplitude of 370 km pointing in the 
direction I = 212°, b = 35°, which is not much different 
from that reported by Faber & Burstein (1988) {v = 360 
km , I = 199°, 6 = 0°) and is greater in amplitude 
than that obtained by Han & Mould (1990) {v = 240 
km s-i , ? = 205°, b = 11°). 

In our minimizing scheme Mark II data favour fio-values 
close to 1, whereas Mark HI subsamples tend to give lower 
rio-values (the whole Mark HI and the Mark HI spirals give 
rio = 0.4±0.3 and = 0.6±0.2, respectively). Although 
the exact estimate of VIq may be seriously affected by our 
specific modeling of the density field in which the detailed 
features of density fiuctuations are neglected, the tendency 
towards lower rip-values is meaningful and agrees qualita- 
tively with the coming down of the estimates of the param- 
eter /? resulting from the comparison between the (Mark 
II and Mark III) observed velocity field and that predicted 
by the IRAS galaxy densities (/3 ~0.4-0.9, with the lowest 
values favored by the most recent works; see, e.g., Sigad 
et al., 1998). Notably, a comparison between the IRAS 
galaxy density fields and the mass density fields recon- 
structed from the application of the POTENT method to 
Mark II (Dekel et al., 1993) and Mark III (Sigad et al., 
1998) data shows a similar tendency, with systematically 
greater values of f3 {f3 ~1.3 and ^0.9, as reported in the re- 
spective papers). Reasons for the systematic differences in 
/?- values obtained from comparing fields at density and ve- 
locity levels have been discussed by Willick et al. (1997b) 
and Willick & Strauss (1998). Direct comparison of den- 
sities is a local procedure, whereas comparison of veloci- 
ties is non-local, because the velocity field is sensitive to 
the mass distribution in a large volume, even outside the 
sampled volume (e.g., the Shapley concentration). This 
non-local character is also inherent to our approach and 
may cause our results to better agree with those resulting 
from velocity-velocity comparisons. 

Also very recent papers, which followed this kind of ap- 
proach, favored fairly low values of /3, by using new sam- 
ples of velocity field tracers. Da Costa et al. (1998) found 
P ~0.6 from an extensive I-band TF sample of spiral galax- 
ies, the SFI catalog (e.g., da Costa et al., 1996). Riess et 
al. (1997) obtained /3 ~0.4 from a sample of Type la su- 
pernovae. 

In conclusion, our best estimate of converges towards 
the (inelegant) value of ~0.5, which roughly corresponds 
to the average of many estimates of flo resulting from the 
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analyses of large scale structure and cosmic flows. Studies 
based on non-linear dynamics within galaxies, groups, and 
clusters (on scales of ^ l — 10h~^ Mpc) yield low values of 
Qo ^0.2-0.3 (see the review by Dekel, Burstein & White, 
1997). 

4 Inverting the redshift-distance relation: 

THE triple-valued REGIONS. 

We use three basic models derived in the previous sec- 
tion, i.e. the modified cluster dipole model and the four- 
attractor models relative to the Mark II subset and the 
Mark III spirals, for providing numerical and explicit ex- 
pressions of the recession velocity cz (expressed in the LG 
frame) as a non-linear function of distance r. In gen- 
eral, this relation, coupled with the heliocentric redshift 
tabulated in Garcia et al. (1993) and transformed by 
us in the LG frame according to the relation czlg ~ 
czq — 79 cos I cos b + 296 sin I cos 6 — 36 sin b (with cz in 
km ) (Yahil, Sandage & Tamman, 1977), gives the 
galaxy distance. The galaxy groups arc given the me- 
dian values of the coordinates and redshifts of the group 
members, selected according to the final catalog by Gar- 
cia (1993). As already said in §3.2.1, because of the Local 
Anomaly problem, for the four-attractor model we do not 
invert the redshift-distance relation for the few very nearby 
galaxies (cz < 700 km s"-'^ ), for which we simply take the 
tabulated distances (with Hq = 75 km s~^Mpc~^). 

It is known that in the vicinity of prominent over- 
densities the redshift-distance relation can become non- 
monotonic, such that there are three difi'erent distances 
corresponding to a given redshift (see the S-shaped curve 
of Fig. 14a in the following section). This makes ambigu- 
ous the distance assignments for the galaxies which fall in 
these triple- valued regions. 

We treat as triple- valued zone objects also the field 
galaxies and the groups which would fall in these zones, 
if their redshifts were modified by ±100 km s~^ and by a 
value equal to the group velocity dispersion, respectively. 
We calculate the velocity dispersions of groups using the 
robust scale estimator defined in terms of the "median ab- 
solute deviation" (Beers, Flynn & Gcbhardt, 1990). In our 
cases all triple- valued zone field galaxies are spirals. 

In practice, if we take the modified cluster dipole model, 
this problem concerns only 4 groups and 14 field galaxies 
in the Virgo region. If we take the attractor model rel- 
ative to the Mark II subsample, this concerns 10 groups 
and 28 field galaxies located in the Virgo or GA regions. 
For the attractor model fitting the Mark III spirals, the 
redshift-distance relation remains monotonic in the vicin- 
ity of GA and there are 10 field galaxies (and the Virgo 
cluster) located in the triple- values regions of Virgo or PP. 

In order to solve this problem, Tully & Shaya (1984) 
simply relied on some DIs, optical appearance, and con- 
nection with some Virgo clouds to guess the distance of the 
Virgo triple-valued zone galaxies. Yahil et al. (1991) used 
a distance-averaging procedure, which however, tends to 
place all objects in the middle branch of the triple-valued 



zone; on the other hand, iteration schemes for translating 
from redshift space to distance space tend to place objects 
outside the middle branch (Yahil et al., 1991; Hudson, 
1993a). Sigad et al. (1998) used a statistical approach 
which takes into account prior information on the peculiar 
velocity field. 

We use the following precepts to solve the problem of 
the triple- valued zones. If a galaxy which falls in these 
zones is included in Mark III, among the three possible 
distances we choose that which is closer to the value given 
by Mark III. If a spiral galaxy is not included in Mark III, 
but has known apparent magnitude Bt and known max- 
imum rotational velocity Vm, we can evaluate the three 
values of the absolute magnitude 

MB-51og/i = Br-51ogr-15.62 (8) 

which correspond to the three choices of distance r (ex- 
pressed in km s^^ ). Then we compare the values of Vm 
and the three values of Mb — 5 log h with the TF relation 
calibrated on a sample of galaxies lying at similar redshift, 
but outside the triple-valued zone. More precisely, for the 
calibration of the TF relation expressed as the linear rela- 
tion 

Mb -51og/i = a- HogKn (9) 

(with Vm in km ), we consider a sample consisting of 
all the galaxies which lie in a distance shell centered on the 
attractor, with the shell having the same width in distance 
as the distance range in which the redshift-distance rela- 
tion is not monotonic. We find the TF relation by means 
of a regression analysis in which we calculate the ordinary 
least-squares bisector line, excluding iteratively the points 
which depart by 3cr from the best regression line. After 
having derived the TF relation, for a galaxy located in the 
triple-valued zone we choose the distance which implies a 
value for the absolute magnitude closer to the TF relation. 

We wish to point out that in this procedure we do not 
want to find the "true" TF relation to be used as an unbi- 
ased DI with low magnitude scatter; we simply use the TF 
as a criterium for choosing the distance of a triple- valued 
zone object, which is found only in restricted volumes close 
to prominent gravitational sources. In these circumstances 
systematic biases act in the same way on all galaxies and 
therefore should not affect our distance choices. 

Table 7 reports the parameters of the TF relations (slope 
b, zero-point o, rms scatter a (in mag)) we used for vari- 
ous regions and the three aforementioned peculiar velocity 
models. The; slopes are in substantial agreement with the 
typical values (6 ~5-6) reported in the literature for blue 
TF relations (e.g., Tully & Fisher, 1977; Bottinelfi et al., 
1983; Garcia et al., 1993). As is expected, the rms scatters 
a about the mean relations tend to be larger than the typ- 
ical value of fT ~ 0.5 — 0.6 mag for blue TF relations (Tully 
& Fisher, 1977), because of errors with which velocity field 
models reproduce the distance, and, especially, because of 
the inclusion in the sample of objects with inaccurate pho- 
tometric and line width data (e.g., galaxies with uncertain 
inclination corrections of the 21 cm line widths). 
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Fig. 13. — Plots of Mb — 5 log h (where Mb is the blue absolute magnitude) versus log VWi (where Vm is the maximum rotation 
velocity expressed in km s~^ ) for the galaxies located in the Virgo redshift region (e), (with galaxy distances predicted by the 
modified chistcr dipolc model), for the galaxies located in the Virgo (c) and GA (d) redshift regions (with galaxy distances predicted 
by the model fitted to the Mark 11 subset), for the galaxies located in the Virgo (a) and Perseus-Pisces (b) redshift regions (with 
galaxy distances predicted by the model fitted to the Mark III spirals) . 



The substantial reliability of the velocity field models 
considered is supported by the fact that they in general 
imply reasonable TF relations. 

Fig. 13 shows the Mb — 51og/i versus logKn plots for 
galaxies located in various redshift intervals corresponding 
to the the triple- valued zones, with distances predicted by 



different velocity models. 

If we had also field ellipticals (which we do not have) 
lying in the triple- valued zones, we could have followed a 
similar procedure, using the modified Faber-Jackson rela- 
tion instead of the TF relation. Remarkably, for the few 
triple- valued zone objects included in Mark III, our 



TABLE 7 

The parameters of the Tully-Fisher relation. 



Model 


Region 


N 


a 


b 


aimag) 


Multi-attractor (Mark III*) 


Virgo 


641 


-6.70 ±0.35 


6.27±0.16 


0.83 


Multi-attractor (Mark 11* ) 


Virgo 


665 


-7.34 ±0.32 


5.93±0.25 


0.79 


Modified cluster dipole model 


Virgo 


743 


-7.28 ±0.40 


5.87±0.20 


0.94 


Multi-attractor (Mark III*) 


PP 


580 


-11.66±0.46 


4.11±0.21 


0.63 


Multi-attractor (Mark II*) 


GA 


957 


-13.15±0.38 


3.39±0.18 


0.70 
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method for choosing the best distance yields the dis- 
tance closer to the value tabulated in Mark III. 

5 Comparing the results relative to different 

velocity field models 

In Fig. 14 we show the predicted behaviors of the re- 
cession velocity and peculiar velocity against the distance 
from the attractor center, along the lines of sight of Virgo, 
GA, and PP. We also show the behavior of the peculiar 
velocity along the line of sight of SH. Fig. 14 and the pre- 
vious Figs. 1, 2, 4, 5, well illustrate the main differences 



between the models. 

There are outstanding features common to the velocity 
field models considered, such as the presence of the attrac- 
tors Virgo, GA, PP and SH; but their prominence appre- 
ciably differs among the models (see the above-mentioned 
plots and §3.2.2). 

Notwithstanding the simple geometry adopted for the 
attractors, our multi-attractor model fitted on Mark III 
data delineates a complex velocity field, which well re- 
sembles the Mark III-POTENT velocity field presented 
by Dekel (1994, 1997), except for the Coma supercluster 
region. 
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Fig. 14. — These plots present the behaviors of the recession velocity cz and peculiar velocity Vp (in the CMB frame), as predicted 
by some models of the peculiar velocity field, against the distance r (in km s^^ ) from the attractor center, along the lines of sight 
of Virgo {I = 284°, 6 = 74°), Great Attractor (/ = 309°, b = 18°), and Perseus-Pisces (/ = 120°, b = -30°). We also show the 
behavior of the peculiar velocity along the line of sight of the Shapley concentration (Z = 308°, 6 = 4°). The dashed, dotted, and 
solid linos refer to the predictions of the modified cluster dipole model and multi-attractor models fitted to the Mark II subset and 
the whole Mark III, respectively. The Hubble relation C2 = r is also indicated. 
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Compared to Mark II, Mark III gives partieularly pro- 
nounced PP and SH attractors and a less dominant role 
for GA in shaping the velocity field. SH appears to ac- 
count for nearly half the LG peculiar motion. Compared 
to the Mark III spirals, the velocity field of the whole Mark 
III, where E/SO galaxies are included, is characterized by 
a greater region which is dynamically dominated by GA, 
but the value of /dof and the errors on the model pa- 
rameters are greater. In Mark III velocity maps there is 
no clear evidence of an infall pattern around PP, which 
rather acts in slowing down motions towards GA in the 
region comprised between LG and PP. On the other hand, 
large streaming flows in the general direction of the CMB 
apex are particularly remarkable. 
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Fig. 15. — We show the compaxisons between the uncorrected 
galaxy distance r = cz (expressed in km and in the CMB 
frame) and the three corrected galaxy distances (hereafter ri, 
r2, rz), as respectively predicted by the modified cluster dipole 
model, the multi-attractor models relative to the Mark II sub- 
sample and the Mark III spirals. 

As regards the cluster dipole model, the recon- 
structed cluster density field (smoothed with a 1500 



km s ^ gaussian window) appears to be in satisfactory 
agreement with the Mark III POTENT density field in 
the nearby region [cz < 5000 km s~^ ) of the supergalac- 
tic plane, where the GA and PP superclusters appear as 
prominent density peaks separated by an extended under- 
dense region (see Plionis et al., 1996). There is less agree- 
ment in farther regions (just where peculiar velocity data 
are less reliable) so much that the Coma supercluster is re- 
placed in Mark III POTENT density maps by an overdense 
region located in a displaced position. The corresponding 
velocity field maps bear less similarity than the density 
field maps, because the modified cluster dipole model is 
characterized by a well-defined infall pattern around PP, 
small large scale flows towards GA and SH, and large Virgo 
infall velocities in the local region. This large infall , which 
is at variance with most recent relevant estimates (see end 
of §3.1.2), enhances the flows towards GA in the LS re- 
gion and makes the infall pattern around Coma almost 
disappearing. Thus, this large Virgoccntric infall makes 
the modified cluster dipole model less dissimilar from the 
Mark III multi-attractor model in the LS region than the 
original cluster dipole model. 

Notably, the picture which is emerging from the appli- 
cation of the POTENT machinery to the new data of an 
ongoing peculiar velocity survey, the I-band TF distance 
survey of about 2000 spiral galaxies in the field (the SFI 
survey) and in the direction of 24 clusters (the SCI survey) 
(Giovanclli et al., 1997 a,b), bears some features which 
arc seen in the modified cluster dipole model, such as a 
clear infall pattern around PP and the absence of coher- 
ent streaming flows on large scale (da Costa et al., 1996). 

Differences in peculiar velocity models give rise to dif- 
ferences in the predicted galaxy distances. Fig. 15 shows 
the comparisons between the uncorrected distance r = cz 
and the three corrected distances (hereafter ri, r2, rz) we 
derive for all fleld galaxies and groups of our sample, as re- 
spectively predicted by the modified cluster dipole model, 
the multi-attractor models relative to the Mark II sub- 
sample and the Mark III spirals. The versus ri plot 
reveals a poor correlation especially at low values, where 
ra > ri in many cases (mostly because of the large Virgo 
infall velocity which characterizes the first model) . The r-i 
versus ri plot reveals a poor correlation also at high val- 
ues, mostly because of strong differences in the infall to- 
wards GA. The ra versus r2 plot shows pronounced devia- 
tions from the one-to-one relation at low values (r2 < 2000 
km s~-^ ), where r2 < in many cases (mostly because 
of differences in the predicted Virgocentric infall) and at 
high values (r2 > 4000 km s~^ ), where r2 > rs for many 
objects (mostly because of the smaller back-side infall to- 
wards GA given by Mark III with respect to Mark II). 
But this plot displays, on average, a smaller scatter than 
the two aforementioned diagrams, especially at low and in- 
termediate values. The comparison between corrected and 
uncorrected distances shows that in general corrections are 
more important at low values, as expected; however, the 
r2 versus r plot displays a large scatter and marked sys- 
tematic effects also at great values. 
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6 Conclusions 

In this paper wc provide homogeneous estimates of dis- 
tances for the individual 3689 galaxies and the 485 groups 
(which contain a total of 2703 galaxies) of a large all-sky 
optical galaxy sample (Garcia et al., 1993; Garcia, 1993), 
which is limited to a depth of 5500 km s~^ and to the 
magnitude limit of completeness of Bt = 14 mag (for 
|6| > 20°). This is the widest, complete, all-sky optical 
galaxy sample for which refined galaxy distances are avail- 
able. All our distance estimates are available on request. 

We recover the distances of our objects by correcting 
redshift-distances for peculiar motions through the appli- 
cation of some peculiar velocity field models. 

We invert the distance-redshift relations relative to dif- 
ferent velocity models and solve the problem of the triple- 
valued zones of this relation by using blue TuUy-Fisher 
relations calibrated on suitably defined samples of objects 
having distances predicted by peculiar velocity models. 

In our work we avoid taking inhomogeneous redshift- 
independent distances obtained from various classical DIs. 
Homogeneous redshift-independent distances are available 
for limited samples of galaxies and it is notoriously diffi- 
cult to combine together the results coming from different 
DIs to achieve an uniform scale of distance. Moreover, 
we do not attempt to correct redshift-distances by using 
the peculiar velocity field derived from the positions and 
redshift of the galaxy sample itself, because this involves 
delicate strategies to recover the galaxy distribution in the 
unsampled zone of avoidance and because our galaxy sam- 
ple does not include all the relevant gravitational sources 
for local peculiar motions. 

We regard the peculiar velocity models considered as 
representative of current common views on the kinemat- 
ics of cosmic flows in the nearby universe. However, re- 
finements of the proper calibration of the Mark HI TuUy- 
Fisher relations, on which our velocity field models de- 
pend, are possible. 

This point was raised by Davis, Nusser & Willick (1996). 
The authors performed a mode-by-mode comparison of 
the peculiar velocity fields expressed in sets of indepen- 
dent basis functions, that were fitted in redshift space to 
the inverse Tully Fisher data from Mark HI and to the 
predictions based on the 1.2 Jy IRAS density field (Fisher 
et al., 1995). Their best-fitting value of /3 ^0.6 yielded 
no acceptable agreement between the two velocity fields 
within a region of 6000 km s^^ radius. On the other 
hand, using a new method, called VELMOD, for maximiz- 
ing the likelihood of the TuUy-Fisher observables, Willick 
et al. (1997b) found that the IRAS-predicted velocity field 
provided a satisfactory fit (with (3 ~0.5) to the Mark III 
Tully-Fisher data, within 3000 km s~^ and on a small 
smoothing scale of 300 km s~^ , if one added an external 
quadrupole, which was essentially expected from the way 
in which the density field was smoothed. The discrepancy 
between these two approaches in comparing velocity fields 
might be a result of systematic errors incurred in match- 
ing the Mark III data subsets, an eflfect to which the latter 



approach is insensitive. 

Very recently, Willick & Strauss (1998) applied an im- 
plemented version of the VELMOD method to an ex- 
panded sample which comprises nearly all Mark III non- 
cluster spirals to 7500 km . The authors confirmed 
that the IRAS-predicted velocity field, with quadrupole, 
was a good fit to the Tully-Fisher data of that sample, for 
P ~0.5 and smoothing scales of 300 or 500 km s~^ . But 
they recognized that the VELMOD Tully-Fisher calibra- 
tion differed significantly from the Mark HI TuUy-Fisher 
calibration for the relatively distant {cz >3000 km s^^ ) 
galaxies of the WCF subsample, which covers particularly 
the region of the Perseus- Pisces supercluster. For these 
objects, the VELMOD TuUy-Fisher relations, which, how- 
ever, rely on the accuracy of the IRAS-predicted peculiar 
velocities, yield distances ^8% shorter than the Mark HI 
calibrations. The main effect is a reduction of the infall of 
the Perseus-Pisces region and, hence, of the bulk flow. 

Ongoing observations aimed at giving reliable North- 
South homogenization (e.g., Strauss, 1997) will cast light 
on doubts concerning the validity of the Mark III Tully- 
Fisher calibration. 

Certainly, major developments in peculiar velocity stud- 
ies will arise by the turn of this century from analyses of 
peculiar velocity catalogs with data of superior quality (in 
terms of sky coverage, accuracy, and homogeneity), which 
can come out from peculiar velocity surveys in progress, 
based on the surface brightness fluctuation (Tonry et al., 
1997), Type la supernovae (e.g., Riess et al., 1997), I-band 
Tully-Fisher (e.g., the SCI and SFI samples by Giovanelli 
et al., 1997a, b), and elliptical fundamental plane (e.g., 
J0rgensen, Franx & Kjaergaard, 1996; Wegner et al., 1996; 
Saglia et al., 1997) methods. New corrections to galaxy 
distances for peculiar motions can be easily implemented 
with the aid of some techniques used in this paper. 

In any case, the use of different models of the peculiar 
velocity field allows us to check to what extent differences 
in the current views on the cosmic flows affect the recov- 
ering of galaxy distances in the nearby universe. We find 
that differences among distance estimates appear to be 
less pronounced in the ^2000-4000 km s~^ distance range 
than in farther or nearer regions (see end of §5). 

These differences can affect the optical luminosity func- 
tion of a galaxy sample restricted to a fairly narrow solid 
angle. For instance, the addition of a Virgo infall alone of 
a few hundred km/s can brighten the characteristic mag- 
nitude M* of the Schechter-type luminosity function by 
a few tenths of magnitude , for a shallow galaxy sam- 
ple (e.g., Efstathiou, Ellis & Peterson, 1988). However, it 
can be proved that these differences have a small effect on 
the optical galaxy luminosity function of galaxy samples 
which, like our sample, cover a very large solid angle (see 
Marinoni et al. 1998b for a detailed discussion). 

In general, peculiar velocity gradients are not strong so 
that all galaxies of a given region share similar peculiar mo- 
tions. But in many regions of the nearby universe positive 
(or negative) peculiar radial velocities of a few hundred 
km/s are present; in these regions our corrections to galaxy 
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distances draw the objects which are nearby in angular po- 
sition closer to each other (or farther from each other) by 
a few Mpc with respect to their uncorrected spatial loca- 
tions. Remarkably, this effect has a minor impact on the 
local galaxy density on large scales (i.e., greater than a few 
Mpc), whilst it affects seriously the local galaxy density on 
small scales (i. e., roughly smaller than 1 Mpc) (sec, e.g., 
Marinoni et al., 1998a for a preliminary account). More- 
over, particularly in the above-mentioned nearest and far- 
thest regions of the volume considered, the latter density 
parameter will be particularly sensitive to differences be- 
tween the various sets of corrected distances. 

The local galaxy density on small scales is an impor- 
tant parameter to be used in statistical studies of environ- 
mental effects on the properties of nearby galaxies, since 
it provides a wcU-dcfincd characterization of galactic en- 
vironment irrespective of membership in galaxy systems. 



Remarkably, much of the observed evolution of the prop- 
erties and populations of galaxies which has occurred dur- 
ing recent epochs {z < 1) can be ascribed to interaction of 
galaxies and their local surroundings. 
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